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Using the covariant spectator theory (CST), we present the results of a valence quark-diquark 
model calculation of the nucleon structure function f(x) measured in unpolarized deep inelastic 
scattering (DIS), and the structure functions g\{x) and g2(x) measured in DIS using polarized 
beams and targets. Parameters of the wave functions are adjusted to fit all the data. The fit fixes 
both the shape of the wave functions and the relative strength of each component. Two solutions 
are found that fit f(x) and gi(x), but only one of these gives a good description of g%{x). This 
fit requires the nucleon CST wave functions contain a large D-wave component (about 35%) and a 
small P-wave component (about 0.6%). The significance of these results is discussed. 



I. INTRODUCTION 

The first measurements of deep inelastic scattering 
(DIS) from polarized protons produced a surprising result 
[TJ [2] which became known as the proton spin crisis [21 2] : 
it turned out that the structure function g\{x) (where x 
is the Bjorken scaling variable, defined below) gave a re- 
sult much smaller that expected. At large Q 2 (where 
q 2 = —Q 2 is the square of the four momentum trans- 
ferred by the scattered lepton), recent measurements at 
a number of experimental facilities 5, 6 give 

r?= f dxg^ix) = 0.128 ±0.013 (1.1) 
Jo 

while theoretical calculations based on the naive assump- 
tion that the nucleon is made of quarks in a pure relative 
S-state give much larger values (Jaffe and Manohar [4] 
give 0.194, and our model gives 0.278, as discussed in 
Sec. |III| below. Note that the experimental value of 0.128 
is remarkably close to the older value of 0.126 [H[2] cited 
by Jaffe and Manohar.) 

The controversy was sharpened by Ji [7J [5] who in- 
troduced a gauge invariant decomposition of the nucleon 
spin into spin and angular momentum components. This 
spin sum rule can be written [SI 110] 

\ = \V + L q + J G (1.2) 

where £ is the contribution from the quark spins, L q the 
contribution from quark orbital momenta, and Jg the 
total gluon contribution. Some experimental estimates 
suggest that £ ~ 0.3 requiring that most of the explana- 
tion for the proton spin come from the other contribu- 
tions, but there is evidence that the gluon contributions 
are small and it is unclear how to interpret this spin sum 
rule [HI- 



With this spin puzzle as background, we decided to 
see what our covariant constituent quark model, based 
on the covariant spectator theory (CST) [T2HT4] . would 
predict for the DIS structure functions. This model was 
originally developed to describe the nucleon form factors 
[15] . and has since been used to describe many other 
electromagnetic transitions between baryonic states, in- 
cluding the N ->• A PH [UJ, the A form factors [18]- 
[20], N -> N* (1440) [2T], N -> AT* (1535) [22], and 
N -> A* (1600) [23]. All of these calculations use con- 
stituent valence quarks with form factors of their own 
(initially fixed by the fits to the nucleon form factors), 
and then model the baryon wave functions using a quark- 
diquark model with a few parameters adjusted to fit the 
form factors and transition amplitudes. The diquarks 
have either spin-0 or or spin-1 with four- vector polariza- 
tions in the fixed-axis representation |24j . 

One shortcoming of our model, as it has been applied 
so far, is that we have not yet included a dynamical cal- 
culation of the pion cloud. Recently, constraints on the 
size of the pion cloud were obtained from a study of the 
SU(3) baryon octet magnetic moments .25, 25], and it is 
clear from this study that the pion cloud contributions 
to the nucleon form factors are not negligible. Further- 
more, even without pion cloud effects it is difficult to 
untangle the form factors of the constituent quarks from 
the "body" form factors (which depend only on the wave 
functions of the nucleon). We need a way to determine 
the nucleon wave functions independent of the contribu- 
tions from the pion cloud and the constituent quark form 
factors. 

The study of DIS provides an ideal answer to this prob- 
lem. In CST, the DIS structure functions directly de- 
termine the valence part of the nucleon wave functions, 
giving both their shape and their orbital angular momen- 
tum content. Adjusting model wave functions to fit the 
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FIG. 1. (Color on line) Feynman diagram for the DIS total cros 
section. All of the intermediate quarks are on shell. 



DIS data fixes all of these components. The calculation 
of the DIS structure functions is also of great interest in 
itself. With this model we can address the nucleon spin 
puzzle directly. 

For this reasons we have decided to "start over" and 
let the valence part of nucleon wave functions be com- 
pletely determined by a fit to the valence part of the DIS 
structure functions. Once the wave functions have been 
determined in this way, the low Q 2 pion cloud contribu- 
tions can be calculated and the nucleon form factor data 
can be used to fix the only remaining unknown quanti- 
ties: the constituent quark form factors. This is planed 
for future work. 

The remainder of this paper is divided into five sections 
and three Appendices. In Sec. [TTJ the DIS cross section 
and structure functions are defined, and theoretical re- 
sults for the DIS structure functions are reported (the de- 
tailed calculations of these structure functions are given 
in Appendix [A]). The calculations use nucleon wave func- 
tions defined in the accompaning paper |27j . and sum- 
marized in Sec. Ill Bl This section also discusses how the 
wave functions are related to the structure functions. In 
Sec. lIIII the data is discussed and it is shown how a model 
without angular momentum components fails. Then, in 



Sec. IV we show how either P or D-state components 
can fix gi(x), and a detailed fit to both the unpolarized 
structure function f(x) and the polarized structure func- 
tion g\(x) is given. We find two solutions, but in Sec, [vj 
we show that only one of them gives a good account of 
the smaller transverse polarization function 172 (%)■ Our 
conclusions are presented in Sec. |VI| Some other details 
are discussed in the remaining Appendices. 



II. STRUCTURE FUNCTIONS FOR DIS 

A. Cross section in the CST 

The DIS cross section can be calculated from the imag- 
inary part of the forward handbag diagram shown in 
Fig. [T] The cross section depends on the hadronic tensor 



d A p' d k m q 



(2nf2E s e q 

= 2tt 



2 9/J.V ) Wl + P^Py 



M 2 



-I\(Gi + G 2 ^fj+X2G 2 \ (2.1) 



where P, q are the four-momenta of the nucleon and the 
virtual photon, respectively, 



Ti-- 

I2-- 



M 
1 



(S ■ q) ie^ a p q a j 



(2.2) 



p' 2 is the energy of the on-shell quark in 



the final state, E s = J m 2 + k 2 is the energy of the on- 
shell diquark with mass m s , and Wi,W2,Gi , an d G2 
are the DIS structure functions. Also in Eq. (2.1) s q ,A 



and A are spin projections of the quark, the diquark and 
the nucleon, respectively, while S stands for the nucleon 
four-vector polarization. 



The hadronic current (J^a)/- 1 is 

(J S A \)» = -u(p',s q )3^q)y AX (Pk) 



(2.3) 



with u(p',s q ) the Dirac spinor for the quark, and j^q) 
the elementary quark current with a gauge invariant sub- 
traction 



i/*(ff) = i«(7/*--Jr) 



(2.4) 



and j q the quark charge operator and "Faa the nucleon 
wave function, with the general form 



*AA = A U(P, A) 



(2.5) 



to be specified below. We emphasize that, in this context, 
the choice of the current (2.4 1 is purely phenomenologi- 
cal, but at least it has been shown in one special case 
that the subtraction term dq^/q 2 arises naturally from in- 



teraction currents neglected here. Once the form (2.4) is 
assumed, it is not necessary to explicitly calculate the 
contributions from the subtraction terms ^q^/q 2 because 
they can be reconstructed from the 7 M term, as discussed 
in Appendix [A] The overall factor of 3 multiplying the 
hadronic tensor arises from the contributions of the three 
quarks [THIHZ] ■ Note that we do not average over the spin 
projections of the target nucleon; this leads to the inclu- 
sion of the last two terms in the hadronic tensor that 
arise when the nucleon target is polarized. 
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The polarization of the nucleon is described by the 
four-vector polarization S with the properties 



-1, S-P = 0. (2.6) 
For a nucleon at rest, this polarization vector is 





sin 9 cos < 
sin 9 sin c 
cos 9 



We use the helicity basis to describe the nucleon spin; 
for a nucleon at rest we will choose the spin axis to be in 
the +z direction, (the direction of the three- vector q). A 
nucleon polarized in the direction of S can be written as 
a linear combination of helicity states, so that 



u(P,S m ) = v[ i) m (<t>,9,0)u(P,X) 



where T> is the rotation matrix 



p(%,0,O) = 



and sum over A is implied. Note that 

(er • S) u(P, S m ) = 2m u(P, S m ) . 

where m = ±^ is the spin projection. 
Using the identity 



'e 2 



u{P,S m )u{P,S m ) = k M {P)\ 



1 + 7° 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



where 



A M (P) 



M + f 
2M 



J2<P,s)u(P,s) (2.12) 



is the positive energy projection operator, and summing 
over the spins of the outgoing quark allows the hadronic 
tensor to be expressed as a trace 



W^(q,P) = 3j2 



d 3 p' d 3 k m q x4 



{2n) 2 2E s e q 



S 4 (p' + k-q-P) 



A j v (q)A mq (p')j^q)0 A A M (P)l[l + 7 5 '$) (2.13) 



xtr 



where the operators 0\ arc defined by Eq. ( |2.5[ ) with the 
wave function spin components given explicitly in the 
next subsection. 



B. Wave function of the nucleon 

In this section we summarize our model of the nucleon 
wave function, which is composed of S, P, and D-state 
components. For details see Ref. [27], hereinafter re- 
ferred to as Ref. I. Briefly, the wave function has a quark- 
diquark structure, with the ith quark off-shell (where 
i = {1,2,3}) and the other two noninteracting on-shell 



quarks treated as a diquark system with total four mo- 
mentum ki and fixed (average) mass m s . The total mo- 
mentum of the nucleon is P. 

The CST wave function of the nucleon is the super- 
position of a leading S-state component, with smaller P 
and D-state components 



(2.7) *Ax(P,k)=n s * S AX (P,k) 



+n P ^ x (P,k)+n D ^ x (P,k). (2.14) 

Each component of the wave function is normalized to 
the same value [see Eq. ( |2.19 ) below] , so if the coefficient 
ns of the S-state is fixed by the coefficients np and np 



(2.15) 



Then the square of each coefficient, n\ (where L = 
{S, P, D}), is proportional to the percentage of each com- 
ponent. The size of the coefficients np and np, will be 
fixed by the fits. The construction of these wave func- 
tions is discussed in detail in Ref. I. 

The formulae are first derived under the assumption 
that isospin is an exact symmetry. The formulae are 
then generalized to allow for the u and d quark distri- 
butions to differ, and all fits were done adjusting the u 
and d quark distributions independently. The S and P- 
state components are a sum of terms with spin (and 
isospin 0) and spin 1 (and isospin 1) diquark contribu- 
tions, while only diquarks of spin 1 can contribute to 
the D-state component. The diquarks of spin do not 
interfere with diquarks of spin 1. The individual compo- 
nents are denoted \I> i '™ with L = {S,P,D} the angular 
momentum and n — {0, 1, 2} labeling the state of the di- 
quark (sometimes the spin=isospin, or as in the case of 
the D-state, all three states have diquarks with spin 1). 
These are 



s,o_ J_ ^.0 
S,l i ,i. 



u u(P,A)Vs(P,&) 



*r ,o = ^0°^(p,A)Vp(p,fc) 

1 ^(eX) a ^(P,A)V^p(P,fc) 



*i?= ^ 0° (e* A ) a G afs (k, Cu)Up(P, X)\k\^ D (P, k) 



aa avTo 



D ,1 
AA = 



1 ^1 J 3 * 



e p D * A Up(P,\)k 2 i> D (P,k) 



Kx = ^t^ 1 (A)pD^{P,k)U a {V,\)^ D {P,k) (2.16) 

where cf) 1 (with / = or 1) are the isospin / parts of 
the wave function (discussed in the next subsection), 
ipL(P,k) are the scalar S, P, or D wave functions, e* A 
is the outgoing spin one diquark four-vector polariza- 
tion with spin projection A in the direction of P, e* DA 
is the outgoing four-vector diquark with an internal D- 
wave structure and spin projection A in the direction of 
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P, and the other four-momenta and operators are 
(P ■ k)P 



k- 



7a = 7q ~ 



M 2 
M 2 



C/ Q (P,A) = ^ 7 5 7 ctU (P,A) 

G a0 (k, Cu) = k a C + - l~g af \k ■ C) 

Dp a {P, k) = kpk a - \gpak 2 

~ _ Pf3P a 
9/3 a - 9l3a - —^p- , 



(2.17) 



with G(k, C, v ) the tensor describing the spin-two coupling 
of a diquark with an internal P-wave orbital angular mo- 
mentum structure to the P-wave motion of the third, off- 
shell quark, and D{P 1 k) the spin- two tensor describing 
the a D-wave motion of the off-shell quark. Note that 
P is orthogonal to all of the polarization vectors, and 
P ■ k = P ■ 7 = 0, implying that P a U a (P, A) = 0, and 



P p D Pa = P a D, ia = 
PbG^^ P a G !3a = 0. 



(2.18) 



As discussed in Ref. I, when isospin is conserved, the 
S, P and D-state wave functions are normalized to 



l=e„ / \MP,k)\ 2 = e (-k 2 )\^ P (P,k)f 

Jk Jk 



= e / k^ D {P,k)Y 

Jk 



(2.19) 



where eo is the charge of the dressed quark at Q 2 = 0. 
With this normalization the square of the coefficients no 
and np can be interpreted as the fraction of the total 
wave functions consisting of D and P-state components. 



C. Quark charge operators in flavor space 



the nucleon (see Ref. I). The isospin operators arc nor- 
malized to 

(0°)^° =1 

E(^)VI=iE t -^ t -^ = 1 ( 2 - 22 ) 



The matrix elements of the quark charge operators in 
DIS are now easily evaluated. Noting that DIS involves 
the square of the quark charge, the result is (including 
the overall factor of 3) 



3 (</>0) T jqjq 4>0 



2\0 



SE^^J^H-SHM^) 1 - ( 2 - 23 ) 



2. Broken isospin 

So far this discussion assumes that the u and d dis- 
tributions are identical, with their relative contributions 
being fixed only by isospin invariance. In fact these dis- 
tributions are quite different at both low and high x, and 
we know that the angular momentum distributions of the 
u and d quarks are also quite different. 

Using the results of Ref. I for broken isospin, Eq. (6.4), 
new effective charge operators (including the overall fac- 
tor of 3) are introduced 



fe 2 ) = 



(e 2 )^|(r^o)(| + iT3)(T.^) = i(| + iT 3 ) 
{el)\ 



!(r.£ + )(§ + §7iO(T.£) 
+ W<-){l + \rs){r-C) 



3U ~ 2 T 3) 



(2.24) 



1. Isospin symmetry 



where it is understood that (e 2 ,) and [e 2 )\ multiply u 
quark distributions and (e 2 )\ multiplies d quark distribu- 
tions. In order to simplify the notation, we write these for 
the proton only; the neutron is obtained from the proton 
by charge symmetry, substituting e u o e^. Hence 



If isospin is a good summetry, the quark charge oper- 
ator in flavor space is 



The isospin operators (ft 1 are 



0°=1 

4>e 



(2.20) 



(2.21) 



(e 2 )° = Se 2 u 
{el)l = el 



(el) 1 



2ei 



where e u = 2/3, 



qjl 

1/3. Hence 



[elf Ml = 3eltftf 



(2.25) 



(2.26) 



where ^ is the isospin three- vector of the diquark with 
isospin projection I. These are operators; to convert 
them to isospin states of the off-shell quark they are mul- 
tiplied from the right by the \ t , the isospin 1/2 spinor of 



with xpL anc i ^pL are ^ ne u an( ^ ^ distributions for the 
L = {S, P, D} state components (see Ref. I). 

We are now ready to report the final results for the 
structure functions. 
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D. Results for the structure functions 

Details of the calculations of the structure functions 
are reported in Appendix [SJ All of the results can be 
expressed in terms of the following functions, which de- 
pend on the quark flavor q — {u,d} and, in some cases, 
on L, which is either the letter {S, P, D} or the number 
£ = {0,1,2} : 



- L (x) = I P 2 (z ) k 2L [^( X )} 2 L > 1 only 



d q (x) = / P 2 (z )fc^,(xX(x), 



and for L = {0, 2} only, 



(2.27) 



where 



9i 



92 



P ■ q 

M 3 



:G 2 



we used the shorthand notation z 
of the structure functions), 

f q =n 2 f q s -2n P n s h° q + 

G P = ^e 2 u g p -e 2 g p ] 



z[x) (where z is any 



■f q p 



G\ — 4 o [29 e u g u 



D 



16 e 2 d g ?] 



(2.31) 



and, for convenience, we introduce the coefficients 



asD 



n s n D 



QpD = -3a / 5 Tip TL]J 



(2.32) 



(x) 



„L+2 



4Mx 



(l-^)^(xXM, (2-28) 



where the integral is 



Mm s 
16tt 2 



dx, 



(2.29) 



to describe the strength of the SD and PD interference 
terms. These results are a summary of the detailed calcu- 



(A83I 



lations leading to Eqs. (A42|, (A45l, (A70), (A76), and 



The formulae (2.301 can be separated into separate u 



and d contributions using the general relations 



(2.33) 



and in these expressions, k — |k| is the magnitude of 
the three-momentum of the spectator diquark (distin- 
guished from the four- momentum k only by context), and 
z is the cosine of t he sc attering angle fixed by the DIS 
condition; see Eq. (2.39) below. The function P2(z) is 
the Legendre polynomial P2{z) — \{iz 2 — 1). Since the 
wave functions depend on only one function x of the four 
momenta P and k [defined and discussed below in Eqs. 
(2.35) and (2.37)], we use the notation ip(P,k) = ip(x) 



for convenience. The physical interpretation of these ex- 
pressions will be discussed in Sec. |IIE| 

In terms of these structure functions, the results for the 
DIS observables for the proton (with the neutron results 
obtained by the substitution e u 4-> e d ) are: 

vWl =2MxWf = x 2el f u + xe 2 f d = xf p (x) 



p 2 1 



f 4 f p 
Jd— 3 Jd 



I 2 2 
l P3 



2 2 r p 



9 a SD [e u d u + 2e d d d ] 



+ la PD [e 2 u h 2 u + 2e 2 d h 2 } 



n 2 P G P 



-^npns 
| clpd 



-n%G p D + l a S D [e 2 , d u + 2e 2 d d d ] 
4elK-h° u ]-e 2 [h d -h d ] 

■ 2 u [hl-h 2 u ]+2e 2 [h d ~h 2 ]\ (2.30) 



with i — {1,2}. Limiting the expansions to u and d 
quarks, and ignoring antiquark, gluon, and correction 
terms coming from the QCD evolution, the extracted f q 
distributions were given in Eq. (2.31 ). The results for the 
<7's are 



9i 



9 d i 



9-2 



fa 



1 •/ u 



■ 5 asD d u 



_ 8 2 f P 
9 n P Ju 

2 1.2 , 

9 apo n u + 



8 T7 2 n P 

9 n P 9u 



^n 2 n D 
60 n D 9u 



If _i_ 4 „2 f P 
~3~Jd "r 9 U P Jd 

-^a SD d d + |a PD h 2 d - %n 2 P g p + j^n 2 D g^ 
\a S Dd u - I npns(hl - h° u ) + la PD (hl 



K) 



4 2 P 
-3 n p9u 



^n 2 n D 
40 n D 9u 



at 



= | a S Ddd + ln P n s (h d - h d ) + ja PD (h d - h 2 ) 



2 2 P 

+ 3 n P9d 



4 2 D 
5 n D9d ■ 



(2.34) 



We conclude this section with a discussion of the inter- 
pretation and normalization of the structure functions. 



E. Physical interpretation 

As discussed in Ref. I [27] , the wave functions are cho- 
sen to be simple functions of the covariant variable 



X 



(M-m s ) 2 - {P-kf 
Mm* 



2P-k 

Mm, 



2. (2.35) 



G 



Since the nucleon and the diquark are both on shell, the 
variable P-k, related to the square of the mass (P—k) 2 of 
the off-shell quark, is the only possible variable on which 
the scalar parts of the wave functions can depend, and \ 
is simply a convenient linear function of this variable. 

For DIS studies we choose to work in the rest frame 
of the nucleon (but our results are frame independent). 
In variables natural to the CST, \ depends only on the 
magnitude of the spectator three-momentum, which will 
be scaled by the mass of the nucleon. If r is the mass 
ratio 



m s 
M 



then k = Mk, E s = ME K = M^r 2 + k 2 , and 



X = x(«)= 



2ME K 



2\/l 



(2.36) 



(2.37) 



The structure functions, displayed in Eqs. (2.27) 



(2.29) (discussed in Appendix \M, are integrals over the 



magnitude of the scaled three-momentum k, or alterna- 
tively integrals over x- I n the nucleon rest frame this 
integral has the form 

Md 3 k 

" 5{E s -kcos9-M{l-x))i> 2 {x) 



{2ir) 3 2E s 
M 2 



(2n) 2 
M 2 



l + x)ip 2 { X ) 



2EL, 



<A 2 (x) 



where the d function in the last integral fixes the scatter- 
ing angle z = cos 9 in terms of the momentum k and the 
Bjorken variable x [the CST form of the DIS scattering 
condition, see Eq. (A39)] 



E K -(l- x) 

K 

r\ + 2(r - 1 + x) 
^Vx(x + 4) 



(2.39) 



The requirement that the scattering angle be physical, or 
that |zq| < 1, fixes the lower limit of the k integration at 



(1 



»- 2(1-,)- (2 - 40) 

for x € [0, 1]. The boundary of the region of integration 
over k, |K m in|, is shown in Fig. [2] for selected values of 
the ratio r. When r < 1, the boundary has a cusp at 
x = 1 — r, which deserves further comment. At the cusp 
l K min| = 0; if x > 1 — r the region k < |K m ; n | is excluded 
because zq > 1 while if x < 1 — r the region k < \ K min \ is 
excluded because zq < — 1. In either case the lower limit 
of x becomes 



1 — x ^ (r + x — l) 2 
r r(l — x) 



(2.41) 



2.00 



1.50 - 



a i.oo 



0.50 



0.00 




FIG. 2. (Color on line) The boundary of the k integral, 
a function of x for selected values of r = 1.25 (short-dashed line), 
r = 1 (solid line), r = 0.75 (long dashed line). 



with (for future reference) 
dx 



2^n 



r(l — x) 

Hence the integral over k can be transformed to 



(2.42) 



M 2 



ndn 
2E~ K 



V> 2 (x) 



Mm s 
4(2tt) 2 



d X ^(x) (2-43) 



(2.38) leading to integrals of the form (2.29) 



It is now easy to interpret the physical meaning of 
the structure functions in the context of the CST. The 
formula 



16tt 2 



dxl^ix)} 



(2.44) 



C 



displays the simplest structure function as an average 
over the square of the wave function [ipg] 2 with a variable 
lower limit that depends on x. The wave function can be 
unfolded from this average by differentation 



[Vf(0] 



r(l - x) 16tt 2 dfljx) 
2K n 

(1- 



m s M dx 
16it 2 dfS(x) 



x(2 — x) m s M dx 



if r = 1 . 



(2.45) 



Examination of the experimentally determined dis- 
tributions, discussed in Sec. |III| below, shows that the 
derivative is singular at x = 0, and is nonzero through- 
out the region < x < 1. Therefore, the wave func- 
tion must have a singularity at x = corresponding to 
k = \\r 2 — 1|. Unless r = 1 this singularity is is at a 
finite and non-vanishing value of k, which is unphysical. 
Furthermore, if we choose r < 1, the wave function will 
have another singularity at x = 1 — r, corresponding to 
C = and k = 0. [In Ref. [15] we used r < 1 and a 
wave function that was finite at all k, inevitably giving a 
distribution amplitude with the wrong shape (see Fig. [9] 
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in Ref. [H]).] We conclude that no choice of r will allow 
us to chose a wave function that is not singular at some 
momentum, but choosing r = 1 where the only singular- 
ity is at k = where Dirac wave functions are known to 
be singular, makes sense physically. With this choice the 
integral (2.29) samples the wave functions over the entire 
range of momentum k, with the sample size depending 
on the value of X , CIS shown in Fig. \2\ 

Assuming for the moment that f£(x) x a (l — x) 1 
at large and small x, this shows that the square of the 
S-state wave function, when r = 1, must go as x a ~ 2 (l — 
x) 1+1 at large and small x. Since 



c 



1 



(2.46) 



this behavior requires the square of the wave function to 
go, at large and small £, like 



l^ s (C)l s 



1 



£l-a/2(£ + 0y+a/2 



where (3 is a range parameter. In this way the asymptotic 
behaviors of the wave function can be estimated directly 
from the structure function, at least for cases where the 
integrand does not depend on the cosine of the scattering 
angle, zq. (Recall that the argument of the wave function 
itself does not depend on z .) We will use this insight in 
Sec.|lV]to guide our construction of models for the wave 
functions. 

We co nclud e this subsection by noting that the ex- 
pression (2.44) display a certain symmetry in which 
is most easily discussed if we introduce y = 1 — x and 
define 



(2.48) 



Then, from the form of £ (written here for r = 1, but 
easily generalized to r ^ 1), 



f"(v) = ft 



(2.49) 



This symmetry allows us to extend the definition of f(y) 
from the interval < y < 1 to the interval < y < oo, 
and will be used in the next section. 



F. Normalization 



The representation (2.44) also allows a simple expres- 
sion for the first moment of fg(x), or fg(y) 



dxf%{x)-- 



Mm* 



Mm, 



dx / d X [^(x)Y 
xdx^ s JOf 



16-7T 2 J dx q 

dyJIiv) 



(2.50) 



This can be compared with the CST normalization in- 
tegral for tf>g . Recalling that the quark charge, at Q 2 = 
is dressed to eg (in units of the charge at high Q 2 ), the 
CST normalization integral is 



d 3 k 



#f(x)] 



e°Mm s 

q a 



8tt 2 



Kdx[^(x)f 



(2.51) 



Note that both (2.50) and (2.51) are integrals over the 



wave function, with \ a function of k and £ a function of 
x. We can transform these two integrals into the same 
form if we set C — Xi which defines a mapping between x 
and k 



r 

T 



r 

1 - x 



1 - x 
r 



(2.52) 



(2.47) Setting r = 1 from here on, the normalization integral 



(2.51) is transformed into 



1 = 



eg Mm s 
16tt 2 

e° Mm s 
16n 2 



X 

1 -x 



dxf^C)] 2 . 



(2.53) 



The normalization integral, apart from the quark charge 
renormalization e°, differs from the first moment of fq (x) 
in the weight function, which now includes the additional 
factor of x/(l — x). 



As an alternative to (2.53), write the integral in terms 



of y, use the symmetry property of £ to transform the sec- 
ond term into an integral over z = 1/y, and then replace 
the integration variable z by y, allowing the normaliza- 
tion integral to be written 



1 = 



16tt 2 
eg Mm s 
16tt 2 



(l-y)dy§-[^(0] 



xdx 



dy 



xdx 



dx 

dfq(x) 
dx 



dxfg(x) 



(2.54) 



in agreement with the results presented in Ref. [TS] . The 
forms (2.53) and (2.54) are equivalent, alternative forms 



of the normalization condition. 

In this model, f q is the valence quark distribution, and 
hence we require its first moment to be unity 



1=1 dxf q (x) 
Jo 



(2.55) 



(where, by our convention, the factor of 2 that accompa- 
nies the first moment of the u quark distribution in the 
proton is contained in the formulae (2.30) and not in the 
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moment). This will set the scale of the wave function. 



The CST normalization condition (2.51 1 will then deter- 
mine the renormalization of the quark charge at Q 2 = 0. 
This is a different procedure than we used in Ref. [T5] . 
and more in keeping with QCD, which fixes the quark 
charges at Q 2 —> oo. 



III. FIRST OBSERVATIONS 

In this section we first present the data for the unpo- 
larized structure functions f q and the polarized g\, and 
then discuss the proton and neutron spin puzzles. The 
detailed fits will be discussed in Sec. IIVI 



G. Where is the glue? 



Data 



Gluons are known to make a substantial contribution 
to the nucleon momentum. This shows up in momentum 
sum rule. Using our normalization, the proton momen- 
tum sum rule is 



dxxf u (x)+ / dxxfd{x) + N g (2.56) 







where N„ ~ 0.5 is the contribution from gluons. This 
sum rule cannot be derived within the context of the CST 
model description of DIS scattering. Instead, we have 
the charge normalization condition (2.54). If the u and 



d distributions are identical (for purposes of discussion), 
and N q = 0.5, these two sum rules are 



0.167= / dxxf°{x) 
Jo 

rl 



dxf*{x) 



(2.57) 



In CST, the first sum rule is fixed phenomenologically 
(by fitting the theoretical f q to experiment) and then 
the dressed charge e q is determined from the second. 

To illustrate the idea, choose the oversimplified distri- 
bution 



f q s (y) = 6(y-y ) + S 



Vo 



(2.58) 



where f q {y) satisfies the symmetry condition (2.491 (as 
it must), and j/o is a parameter. Then, the normaliza- 
tion condition (2.55) is automatically satisfied, and the 
conditions (2.57) become 



xdxf^x) = l-yo = 0.167 



i-S( 1 + 4 



(2.59) 



These equations give e q = 0.41, not too far from the 
values obtained in the following sections. 

It remains to be shown in detail how the gluon contri- 
butions give rise to the modification of the quark charge, 
and to generalize this discussion to show how the angu- 
lar momentum contributions of constituent quarks, eval- 
uated in this paper, can be compared to the angular mo- 
mentum contributions from bare quarks plus gluons that 
would be obtained from a light-front model. This is be- 
yond the scope of this paper and a subject for future 
work. 



The individual u and d quark distributions f q and g\ 
have been extracted from global fits to data. These fits 
use the QCD evolution equations to relate the data at 
higher Q 2 to phenomenological starting distributions de- 
fined at Q 2 = 1 GeV 2 . In fitting our model for f q and 
g\ to these starting distributions we assume that it is 
appropriate to compare the DIS limit of our model with 
the data at Q 2 = 1 GeV 2 , assuming that the behavior 
at higher Q 2 can be predicted by QCD but not by the 
model. Some investigators have chosen to extrapolate the 
QCD predictions to (much) lower Q 2 and fit their models 
there. We do not do so for two reasons; (i) we are doubt- 
ful that the QCD evolution equations are reliable below 
Q 2 = 1 (they are based on perturbative QCD), and (ii) 
the assumption that our model has reached its asymp- 
totic limit at the scale of the nucleon mass is generous. 
The choice of Q 2 = 1 GeV 2 is somewhat arbitrary, and 
at best a compromise required by trying to fit a round 
peg into a square hole. 

For the fits to f q we use the global fits of Martin, 
Roberts, Stirling and Thorne (MRST02) [30] : 



xf? v (x) 
xfT P (x) 



xu v (x)/2 =0.130 x a31 (l-x) 3 - 50 

x(1 + 3.83^ + 37.65:e) 

xd v (x) = 0.061322 x ' 35 (l - x) 

x(l + 49.05yi + 8.65x), (3.1) 



\4.03 



where we have divided their u quark distribution by 2 
because both of our f q {x) distributions are normalized 
to unity 



dxf^(x) 



1 



(3.2) 



[In order to get this valence quark normalization condi- 
tion accurately, we we rescaled the MRST02 from 0.131 
to 0.130 (for u quarks) and from 0.061 to 0.61332 (for d 
quarks). These rescaled numbers are shown in bold in 
Eq. ( |3.1[ ).] The model does not describe sea quarks, so it 
is appropriate to use the valence quark distributions; the 
description of sea quarks is a subject for future study. 

To represent the data for gf, we use the global fits of 
Leader, Sidorov, and Stamenov (LSS10) 31J. They ex- 
press these distributions as a sum of the leading twist 
component, Aq, and a higher twist component propor- 
tional to h q 



g\{x) = Aq{x) + 



h*(x) 



(3.3) 
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These values are compatible with recent measurements 
reported at a number of experimental facilities [S] , and 
also agree, within errors, with the result for the proton 
reported by Jaffe and Manohar 0]. 

For later use, we record the experimental values of T± 
for the separate u and d distributions obtained from the 
experimental results using the expansion (2.33): 



0.2 



0.4 



0.6 



r " = ! ( 4r i - r i ) = °- 333 ± °-° 39 

rf = | Ur^ - r{\ = -0.355 ± o.oso , (3.9) 

where the errors are estimates obtained by integrating 
Eq. ((3J). 



FIG. 3. (Color on line) The data for h N (x) and the empirical fits 
of Eq. l |3.5| . The dashed lines above and below the solid curves are 
the error estimates of Eq. < |3 - 5 [ ) - Lower points and curves are the 
proton; upper are the neutron. 



At Q 2 = 1, the leading twist contributions determined 
by LSS10 are 

xAu(x) =0.548 x°' 782 (l - a;) 3 ' 335 

x(l - 1.779y/x+ 10.2 x) 
xAd{x) =-0.394 a; - 547 (l - x) 4 ' 056 (l + 6.758 ar) . (3.4) 

However, the higher twist corrections are not negligible 
at Q 2 — 1. Extracted values of h N for the proton and 
neutron are reported in Table III of Ref. [31] . For conve- 
nience these were fit by the smooth functions; 



hl(x) 



0.82x v -' az {l-xf[l±S p (x)] 



hl(x) = 3(1 - x) 4 (0.2 - x) 2 [l ± S n (x)} 
where the errors were approximated by 

0.05 0.04 

(1 - .t) 4 + ~x^~ 



(3.5) 



S p (x) 

6 n (x) 



r 0.7 



0.015 



0.04 



(0.2 



(1 



(3.6) 



The quality of these fits and the errors in the h's are 
shown in Fig. [3] Separating these into u and d quark 
contributions and adding them to (3.4), the total empir- 



ical distributions (adding errors) at Q = 1 become 

xgH P± (x) = xAu{x) + lx{4h v ±{x) - h n ± {x)] 
xgZ P± (x) = xAd(x) + lx[4hl(x) - h p ± {x)] . (3.7) 

As Fig. [3] shows, the error in these functions is consider- 
able. 

Forming the combinations for proton and neutron, and 
integrating over x gives 



r?= / dx 9tp P± ( x ) = 0.128 ±0.013 

dx gll p± (x) = -0.042 ± 0.013 . (3.8) 



1 1 



B. Proton and neutron spin puzzles 

Note that, if the nucleon has no P or D-state com- 
ponents, the polarized spin structure functions g™ are 
uniquely predicted. From Eq. (2.30) we obtain 



9i( x ) —30 fp 15 fn — 27 f u 54 f d 



9l ( x ) —is fp 15 fn — 27 f u 27 f d (3.10) 



Recalling the normalization (3.2), the moments predicted 
by J3T0| are 



r?= f dxg p 1 (x) = i= 0.278 
Jo 



r? = / dxg?{x)=0 
'o 



(3.11) 



These predictions are to be compared with the experi- 
mental values (3.8). As mentioned in the Introduction, 



the inability to correctly predict T\ has been referred to 
as the proton spin puzzle, and we see here that our model 
without orbital angular momentum components cannot 
reproduce the experimental results. From our point of 
view, the neutron spin puzzle is also interesting. 

The extent to which the prediction (3.10) for g\ 
strongly disagrees with the data is shown in the left panel 
of Fig. [4] the prediction for g\ shown in the right panel is 
in better agreement. This figure is included here only as 
an illustration of the spin problem. The data shown are 
only a subset of all the data that is available and were 
measured at a variety of Q 2 , which, because of the effects 
of QCD evolution, should strictly speaking not be placed 
on the same plot. The method we used to fit the data, 
which allows for QCD evolution and includes a larger 
data set, is discussed in the next section. 



IV. FITS TO THE DATA 

The fits to the DIS data will be done in three steps. 
First, the unpolarized structure functions f q (x) will be 
fit using a model with only an S-state component. 
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FIG. 4. (Color on line) Data for xg\ compared to the predictions ( |3.10[ l (solid line) and the LSS10 fits (for 5 = 0) l |3.7| (dashed line). 
Left panel: proton with data from SMC (circles, Ref. [32]), SLAC-E143 (squares, Ref. [33]), and HERMES (triangles, Ref. pg]); right 
panel: neutron with data from SLAC-143 (squares, Ref. [33]), HERMES (small circles at low x, Ref. [34]), SLAC-E154 (medium circles, 
Ref. [35]) i JLab-Hall A (large circles, Ref. |360, and JLab-Kramer (diamonds, Ref. I37| ) . 



After the S-state component has been fixed, we use 
the same parameters for the P and D-state components 
and adjust the strength parameters n p an d np, to get the 
correct values of and Tf given in (3.9). Once the np 
and np, have been chosen, we readjust some of the pa- 
rameters of the wave functions to give a good fit to the 
shapes. This last step confirms that our choices of np 
and np, made in step 2 were acceptable. However, this 
procedure is rather crude, and a fit to all of the parame- 
ters at once would likely alter our conclusions somewhat. 



A. Step 1: fitting the S-state wave functions 

As described above, the S-state u and d quark distribu- 
tions are are fit to the experimental quark distributions 
/„ and fd using Eq. (2.44). We choose a simple form for 



the wave functions, with parameters adjusted to give a 
reasonable fit to the data. 

In our previous work |15j we made the choice 

N s 



(4.1) 



m s [x + Pl][x + fa] ' 
where Ns is a normalization constant, Pi and P2 are 
dimensionless range parameters, and m s is included to 
make the normalization constant dimensionless. This 
form, used in our previous analysis of the nucleon form 
factors, does not do well fitting the DIS data, and, as 
discussed in the Introduction, we will adopt a completely 
different approach in this paper. 

A choice that produces a good fit to the DIS data is 

, S( v = 1 Ps q cos9 Sq + xsin9 Sq 

= {m a N?)- x * a q ix,PAno,ni) (4-2) 



where Njjf is a normalization constant, Ps q is the sin- 
gle dimensionless range parameter, 8s q = as q 7r a mix- 
ing parameter which allows the phase and/or oscillations 
of the wave function to be adjusted as needed, n 0Sq is 
a fractional power needed to give the sharp rise in the 
distributions amplitudes at small x, and n\s q allows for 
adjustment of the large x behavior of the wave func- 
tion. Th e ge neric wave function, $(x, P, 9, no, n%), de- 
fined in (4.2), will also be used to define the P and 
D-state wave functions below, and we use the notation 
®«(X, Pr--) = Ps q , ■■■)■ If a Sq ^ and n 1Sq = 3, 



Eq. (4.2) still goes as 1/x 2 at large x, ensuring that the 
form factors calculated from this wave function will go 
as l/Q 4 at large Q. 

Taking r — 1, and using the correspondence presented 
in Eq. (2.47), and assuming 6s q ^ 0, the leading behavior 
of the distribution amplitudes f q near x — > and x — > 1 
should be of the from x aq (l — a;) 7 ', where 



n s q 



2 4°9 



hi- 



(4.3) 



For the u quark distribution, we will fix n\s u = 3, 
giving a (1 — a;) 3 behavior at large x, so this simple model 
cannot reproduce the fractional power of (1 — a;) 3 5 found 



Ps q 



9sa 



nosq 



nisq 



0.9 
1.25 



0.47T 

l _ 



0.51 
0.49 



3 
3.2 



2.197 
2.279 



0.3545 
0.3940 



TABLE I. Adjustable parameters (in bold) and additional con- 
stants determined by the normalization conditions for the fits 
shown in Fig. [5] 
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FIG. 5. (Color on line) The experimental quark distributions 
2;/u Xp (^) and xf^ p (x) (dotted lines) compared to the fits xf^(x) 
and xf^(x) with the parameters given in Table | l| (so lid lines). All 
distributions are normalized to unity as in Eq. | [3.2[ l . The curves 
that peak at a larger x are the u distributions, the ones that peak 
at smaller x are the d distributions. 




FIG. 6. Comparison of the fitted to the experimental ratio of the 
d and u quark distributions, where q(x) = f q (x). The dotted line 
represents the experimental results. 



in the empirical u quark function (3.1 ). However, in order 



to preserve the interesting d/u quark behavior as x — > 1, 
we allow the power nisd > 3, and adjust it to give a good 
d/u ratio. 

If the wave functions were to reproduce the small x 
behavior exactly, we would require 



nosq 



1 - if* 

2 4 iX 1 



0.42 u quark 
0.41 d quark . 



(4.4) 



The actual parameters that emerge from our fits are not 
far from these estimates. 

The S-state structure functions reduce to 



1 



ee {Clr'F^x) (4.5) 

where, for simplicity, the parameters that enter <E>^ have 
been suppressed, and 



(4.6) 



To complete the determination of fg(x), C q and e q < 1 
are determined by the wave function normalization con- 
dition (2.54) [or (2.51)] and the normalization of the va- 
lence quark distribution (3.2). The wave function nor- 
malization gives Cq in terms of e q 



C; 



dxF^(x). 



with e q given by the ratio 



1 

N 



dx f (x 



with 



N 



dxf*(x). 



(4.7) 



(4.8) 



(4.9) 



The fit is shown in Fig. [5] and the parameters are given 
in Table]!] Note that wave function is highly singular (but 
they are normalizable; only if n s q > 0.75 can they not be 
normalized), and even more singular than the estimates 
(4.4). The ratio of the d to u quark distributions, which 
vanish at x —¥ 1, are shown in Fig. [6j We constrained 
nisu = 3, but we cannot get the correct result for the d/u 
ratio unless nisd > 3, and choosing the power nisd — 
3.2 gives the fit shown in Fig. [6j This means that the 
asymptotic nucleon form factors will be dominated by 
the u quark distribution. 



B. Step 2: fixing the P and D-state admixtures 

The strength of the P and D-state admixtures is fixed 
by fitting the moments of the polarized distributions 
(3.8). Before doing this, we discuss our choice of the 



functional form for the P and D-state wave functions, 
and then look qualitatively at the effect of changing np 
and no individually. 



1. Functional form for the wave functions 

We chose the P and D-state wave functions to have the 
general form 

<W = [m s N q D M 2 K 2 { X )]- 1 ^{x^, Mo.ni), (4.10) 

where q — {u,d}, N q are normalization constants, and 
the function K(x) is J us t the scaled three-momentum ex- 
pressed as a function if x 

K(x) = ^lri X ( X + 4) = K, (4.11) 

and the generic functions $^ have the same general 
form as but with different parameters associated 
with each functions [so that, for example, §g(f3, •••)—> 
*(&>„••■)]• 
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FIG. 7. (Color on line) Influence of a D-state on the spin polarized structure functions, xAq(x), as a function of x . In all panels, the 
dotted line is the experimental fit of LSS10, and the theoretical curves have = (solid), 0.2 (long dashed), 0.4 (medium dashed), and 
0.6 (short dashed). Left panel is xAu, middle panel xAd for nrj > 0, and left panel xAd for njj < 0. 



The factors of K 1 multiplying the P-state wave func- 
tion and K ~ 2 multiplying the the D-state wave function 



in Eq. (4.10) were chosen for convenience; they cancel 



similar factors in the spin structures multiplying theses 
scalar wave functions. Our philosophy is to let the short 
range behavior be fixed by the structure functions them- 
selves, which are more singular than might be expected. 

to fix zo, and substituting the g enera l 
;he structure functions that appear in ( 2.27 ) 



Using (2.39 



forms (4.10), t 



and (2.28) consist of terms depending on the squares of 



the wave functions 



d X [^(x)?^(C^- 1 F g L (x) 



dxP 2 (zo)[^(x)} 2 



(4.12) 



and interference terms 



1 f°° 

d q (x) = . / d X P 2 (z )$ S q (x)* q D (x) 



h L Jx) = 



h L q + \x) = 



CfCfrJc 



dXZO ®q{x)®q{x) 



dX^il-zD^ix^ix) 



4.T 



where £ was defined in Eq. (2.41 ), 



C L q = 16vr 2 (A^) 2 , 



(4.13) 



(4.14) 



and the L's appropriate to each equation were specified 
before. The functions h q +1 enter into the expressions for 
<72 only, and will not be needed until the next section. 

The P and D-state wave functions are normalized as 
in ( |2.19 ), and with the ansatz (4.10) this leads to a gen- 
eralization of the condition (4.7) 



C L q = e 



° f dxF q L (x). 

J — OO 



(4.15) 



However, the valence quark normalization condition 
needed to extract e q is now 



1=1 dxf q {x) 

fl 



dx 



n 2 s f q 3 - 2n P n s ti q + n P ff + n< D ft 



D 



(4.16) 



Since the individual wave functions are all normalized to 



the same quantity (e q ) 1 , the constant TV defined in (4.9 ) 
is independent of L 



N = 



dxf q {x) 



dxf?(x), (4.17) 



so that the new result for e q is a generalization of (|4.8) 



p° - n 2 p S 



l S c q 



2n P n s h SP + n 2 P e P + n 2 D ef (4.18) 



where the new coefficients e q are generalizations of (4.8) 



' L q =^J dxft(x) hS q P ^^J Q dxh q (x).(A.l9) 



Next the choice of the coefficients n P and n D is dis- 
cussed. 



2. Effects of the D-state 

First, consider the results of varying the D-state ad- 
mixture with the P-state component identically zero. 



TABLE II. Values of for the four choices of no shown in 
Fig. 



0.2 0.4 0.6 expt 



T1 (n D > 0) 0.667 0.643 0.544 0.367 0.333 
rf(n o > 0) -0.333 -0.293 -0.252 -0.218 -0.355 
rf (n D < 0) -0.333 -0.369 -0.395 -0.404 -0.355 
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Keeping the parameters for the D-state wave function 
identical to those for the S-state (Table]!]), the predictions 



of Eq. (2.34) reduce to (ignoring the small g® terms here, 
but not m the calculations and figures discussed below) 



it 2 tS 2 pD 2 j 
9l - 3/u - n D Ju ~ 9 a SD a u 

d _ 1 f S 8„ j 
9l — -3Jd — g a SD Ud , 



(4.20) 



where, in this case, fq = fq ■ The results for various val- 
ues of no are summarized in Fig. [7] and Table [11} Note 
that, when the D-state parameters are identical to the 
S-state parameters, the normalization of the charge is 
unchanged, i.e. e° = e^. Since the unpolarized distri- 
butions f q {x) are identical to the S-state distributions 
shown in Fig. [5] Fig. [7] does not show f q . 

As the left panel ofFig. [7] and the first line of Table [TT| 
show, a D-state component of about 36% will fit the u 
quark polarization and T" very well. It also shows that 
the sign of np> must be positive in order for the interfer- 
ence term, d u , to bring the predicted g™ down close to the 
observed shape. [A detailed study of Eq. (4.20) reveals 
that a large negative no close to unity is also mathe- 
matically possible, but this case is rejected on physical 
grounds.] Unfortunately, the d quark distribution, fit well 
by n D = 0, is destroyed by increasing \n D \ to 0.6. Since 
the phase of the d quark D-state wave function can be 
chosen independently of the phase of the u quark wave 
function, both possible choices of sign for the d quark dis- 
tribution are shown, but neither will correct the problem. 
The principal source of the difficulty is the relative size of 
the SD interference them, which is four times larger for d 
quarks than for u quarks, and is not compensated by any 
positive term of the form n 2 D f® [the term proportional 



to n D gj[ , neglected in the estimate (4.20), is positive, 



but much to small to have the desired effect]. 

However, as will be shown, a reasonable fit can be ob- 
tained by adding a small P-state contribution and ad- 
justing the parameters of the D-state distributions. 



solution np(u) no(u) 



n P (d) n D (d) 



0.43 
0.08 



0.18 
0.59 



-0.43 
0.08 



-0.18 
-0.59 



TABLE IV. Values of np(q) and no(q) for the two solutions 
studied in this paper. When the sign of n_t(it) 7^ riL(d) it 
means that the phase of the wave function differs from that 
of the S-state. 



to those for the S-state (Table pL the predictions of Eqs. 
(2.31) and (2.34) reduce to (ignoring the small g^ terms 
here, but not ni the calculations and figures discussed 
below) 



/„ = /f - 2n s n P h° 



Si 
9t 



'q 

2 rS 



\n s n P h° u 



^n 2 
9 " 



P fu 

p fd 



til 



(4.21) 



where the SP interference term included in f Q has been 
separated out explicity, and, following Eq. (4.18), the 
renormalized charge reduces to 



= e£ - 2n s n P h° 



SP 



(4.22) 



The results for vari ous values of np are summarized in 
Fig. IS] and Table III Now that e° Q and the unpolarized 



structure functions depend on rip, these quantities are 
also shown in the table and the figure. 

The equations and figure show that the most efficient 
way to reduce the polarized u quark distribution to the 
correct size is to require np > for u quarks, but as the 
phase of the d quark wave function is independent of the 
choice for the u quark, we show both cases. The choice 
of np = —0.4 for the d quark preserves the shape of both 
of the d quark structure functions. We therefore expect 
a solution in the region with np > for u and np < 
for d quarks. 



3. Effects of the P-state 

Now consider the results of varying the P-state admix- 
ture with the D-state component identically zero. Keep- 
ing the parameters for the P-state wave function identical 



4- Fixing ns and no from the moments 

An efficient way to estimate the strength of the P and 
D-state mixing is to adjust the coefficients np and np> to 



TABLE III. Values of r? and e° for the four choices of np shown 
in Fig. H 



\n P \ 





0.2 


0.3 


0.4 


expt 


17 (n P > 0) 


0.667 


0.627 


0.563 


0.448 


0.333 


e° u (n P > 0) 


0.355 


0.269 


0.229 


0.194 




rf (n P > 0) 


-0.333 


-0.313 


-0.280 


-0.223 


-0.355 


e° d (n P > 0) 


0.394 


0.303 


0.261 


0.224 




rf (n P < 0) 


-0.333 


-0.321 


-0.307 


-0.289 


-0.355 


e° d (n P < 0) 


0.394 


0.485 


0.527 


0.564 





case 


A 1 


e ° 


rf 


e° d 


lo 


0.332 


0.187 


-0.356 


0.571 


1a 


0.298 


0.170 


-0.333 


0.567 


2o 


0.334 


0.326 


-0.356 


0.364 


2a 


0.334 


0.326 


-0.294 


0.364 


expt 


0.333±0.039 




-0.355 ±0.08 





TABLE V. Values of T\ and e° for the cases shown in Figs. 
|9"1- |12[ Cases lo and 2o are solutions 1 and 2 with P and D- 
state parameters identical to the S-state parameters of Table 
[I] cases 1a and 2a have the parameters given in Table |VT| 



14 




FIG. 8. (Color on line) Top row: unpolarized structure functions xq = xf q ; second row: spin polarized structure functions, xAq(x), all 
as a function of x. In all panels, the dotted line is the experimental fit of MRST02 or LSS10, and the theoretical curves have \np\ = 
(solid), 0.2 (long dashed), 0.3 (medium dashed), and 0.4 (short dashed). Left panels are xu and xAu, middle panels are xd and xAd for 
rip > 0, and left panels xd and xAd for rip < 0. 





FIG. 9. (Color on line) Distributions for the u quark, as a function of x, for solution 1 with rip = 0.43 and np> = 0.18. Top row: 
unpolarized structure functions xu = xf u ; second row: spin polarized structure functions, xAu(x); left column: same parameters for all 
L (case lrj); right column: parameters readjusted to fit the shapes (case 1^)- 



15 





FIG. 10. (Color on line) Distributions for the d quark, as a function of x, for solution 1 with np = —0.43 and np> = —0.18. See caption 
to Fig. [9] for more details. 




FIG. 11. (Color on line) Distributions for the u quark, as a function of x, for solution 2 with np = 0.08 and np> = 0.59. In this case the 
initial fit was satisfactory, as discussed in the text (case 2q is identical to 2^). 
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fit the experimental moments T\. If the parameters of the 
P and D-state wave functions are set equal to the S-state 
ones, then the theoretical expression for these moments 
gives the following equations for np and nu 



coefficients are 



1 1 



§nl)F, 



)F U — {qCISD 



p 

0.333 

r x = — | + ^np)F d — (fctSD 
+ la PD H d = -0.355 



9 n P 



60 ,l D) Kj u 



±n 2 



j^n 2 D )Gd 



(4.23) 



where we require that the solutions to the first equation 
give positive values of both np and nn, but allow solu- 
tions with either sign for the second equation, and the 



q 


L 




@Lq 


TlOLq 


TllLq 




Ml 





0.7 


0.45tt 


0.45 


3 


2.525 




1 


0.7 


0.47T 


0.45 


3 


1.473 




2 


0.7 


0.47T 


0.45 


3 


1.473 


dx 





1.8 




0.53 


3.2 


5.152 




1 


1.8 


-0.27T 


0.53 


3.2 


5.396 




2 


1.8 


-0.36tt 


0.53 


3.2 


22.153 


U2 





0.9 


0.47T 


0.51 


3 


2.197 




1 


0.9 


0.47T 


0.51 


3 


2.197 




2 


0.9 


0.47T 


0.51 


3 


2.197 


d 2 





0.5 


0.45?r 


0.41 


3.2 


1.082 




1 


0.5 


\- 


0.41 


3.2 


0.234 




2 


0.5 


-0.27tt 


0.41 


3.2 


0.497 



TABLE VI. Parameters for the cases 1.4 (denoted by Mi and 
d\) and 2a (denoted by 112 and (fe); those that were adjusted 
are in bold. The fits are shown in Figs. [9]- |12| 



F q = fdxf^{x) H q = fdxh^\ 
Jo Jo 

G q = f dxd q (x)= I dxg q ia} {x). 



(4.24) 



With these requirements, there are only two solutions 
to these equations, summarized in Table |IV| The values 
of T'J and e q for these two solutions, with the parameters 
for all the P and D-state components identical to the S- 
state, as discussed above, are given in Table [Vj and the 
shapes of the quark distributions are shown in Figs. [9] 
Without further changes in the parameters of the 
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individual components, these solutions do not describe 
the shape of f q and gf, but do give (within round off 
errors) the correct moments T\. 

To complete the fits we move to step 3, adjusting the 
parameters of the wave functions to give a good descrip- 
tion of the shapes. 



C. Step 3: adjusting parameters to fit the shapes 

As a last step, some of the parameters /3l 9 , and 
n 0Lq a re adjusted to fit the shapes of f q and g\ without 
destroying the moments Y\. The results are summarized 
in Table [y] and Figs.|9]-[l2] In all panels of these figures, 
the dotted lines are the experimental fits of MRST02 or 
LSS10 and the shaded band is our estimate of the exper- 
imental errors in g\ obtained from the curves in Fig. [3] 
and Eq. (3.5 1. The solid lines show the full theoretical 
result. 



The panels showing xf q {x) = xq(x) also show 
the individual contributions x{n 2 s f^ + n pf q + npjfq) 
(short dashed line) and —2nsnpxh q (long dashed line) 
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FIG. 13. (Color on line) Data for xg2 compared to the predictions from Eq. | |2.30| for the two solutions. In each panel the total result is 
the long dashed line (solution 1) or solid line (solution 2). The predictions for np = are the dotted line (solution 1) and short dashed 
line (solution 2). Left panel: xgt^, with data from SLAC-E143 (circles, Ref. [33]) and SLAC-E155 (squares, Ref. [38]); Right panel: xg% , 
with data from SLAC-E155 (circles, Ref. [38]), JLab-HallA (squares, Ref. 36 ), and JLab-Kramer (diamonds, Ref. I37| ). 



which add up to the total result xf q . In the panels show- 
ing xg\(x) = xAq(x) the individual results shown are 
those proportional to xf q (the first term on the rhs of 
Eq. ( 2.34[ ) - short dashed line) and those proportional to 
the remaining terms on the rhs of Eq. ( 2.34[ ) , all of which 
would vanish if np = np> = (long dashed line). In the 
left panels of Figs. |9j [lOj and[l2j the parameters of all of 
the component wave functions (S, P and D) are given in 
Table [I] i n the right panels the parameters are given in 
Table rvTl 

Note that the final fits and values of all lie within 
about one standard deviation of the values obtained from 
the best LSS10 fit. The parameters space has not be sys- 
tematically searched: in all cases the values of j3 and no 
are equal for each L and q, and, to the level of accuracy 
achievied here, it was not necessary to refit the distri- 
butions U2- The 9 parameters give great flexibility to 
change the shape of the wave functions, and have been 
used to significant effect, particularly for the d quark D- 
state wave functions. Note that a change in 6 by 7r will 
merely change the sign of the wave function, leaving its 
shape unchanged. This could be used to fix the d quark 
values of np and np> to positive numbers. 

From our study so far, we can only conclude that the 
the two solutions 1 and 2 do equally well in describing 
the structure functions f q and g\. Significant differences 
appear when the predictions for gi are examined. 



PREDICTIONS FOR THE g 2 
DISTRIBUTIONS 



The data for the structure functions gi is incomplete 
and inaccurate but sufficient to allow us to draw some 
interesting conclusions. 



The predictions for g% and g% that follow from the two 
solutions are shown in Fig. 13 The theoretical predic- 
tions are given in Eq. (2.30); note that gi are the only 



structure functions that depend explic itly o n the inter- 
ference terms h q +1 , introduced in Eqs. (2.28) and (4.13). 

The data for g 2 (x) seem to favor the solution z. In 
both panels the total prediction for solution 2 (solid lines) 
gives a very reasonable approximation of the data, even 
though the large error bars do not allow us to draw this 
conclusion too confidently. But it is clear that the total 
predictions for solution 1 (long dashed lines) even fail 
to describe the data qualitatively, particularly for the 
proton. (However, it is interesting that both solutions 
describe the very accurate JLab neutron data equally 
well.) We conclude that solution 2, with its large D- 
state, seems to be a more accurate model. Note that even 
though this model has a small P-state probability of only 
about 0.6%, the P-state does not play a negligible role in 
the description of the data. For example, the left panel 
of Fig. [13] shows that the P-state contribution improves 
the fit to the data for xg^ix) substantially (compare the 
short dashed line with the solid line). 

Another interesting observation is that the theoreti- 



cal formulae, Eq. (2.34), show no sign satisfying or even 



approximating the interesting WW relation |39j 

92{x) = -gi(x) + ( ^giix'). 

J X x 

Numerical studies of the inequalities (for n > 0) 



(5.1) 



dxx n 



n+ 



i9i(x) + g2(x) 



dxx n g 1 (x) (5.2) 



which lead to the WW relation, also show no sign of 
being satisfied by our model. Most striking, perhaps, is 
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that the Burkhardt-Cottingham sum rule 



dx g2{x) = 



(5.3) 



is satisfied only if no = rip = 0, but, because of the 
presence of the structure functions hg +1 singular as x — > 
0, diverges otherwise. Perhaps sea quark contributions, 
neglected here, are essential for a full understanding of 
these realtions. Some of these mysteries may be clarified 
by future study. 



VI. SUMMARY AND CONCLUSIONS 

This paper uses ideas from the Covariant Spectator 
Theory to develop a covariant constituent quark model 
of the nucleon that explains the general size and shape 
of all of the DIS structure functions of both the proton 
and the neutron: f(x) measured in unpolarized scatter- 
ing and gi(x) and gi{x) measured in polarized scatter- 
ing. We find that a good semi-quantitative description is 
possible even thought we treat the larger valence quark 
contributions only, ignoring the smaller sea quark con- 
tributions. The nucleon wave functions include S, P and 
D-state contributions, with the sum of their probabili- 
ties normalized to unity, and the strength of the P and 
D-state contributions proportional to the parameters np 
and no- These parameters, together with those defining 
the shape of the wave functions, are adjusted to fit the 
data for f{x) and gi(x), and the models defined by these 
fits are then used to predict g2{x). 

The major effort required to construct the angular mo- 
mentum components of the wave functions is described 
in detail in the companion paper |27| and summarized 
in Sec. |IIB| In that work we built the general covariant 
forms for the angular dependence of three-body states 
with total angular momentum L — 1,2. Once the wave 
functions have been parameterized, the calculation of the 
DIS cross section is straightforward; the details of the 
evaluation of the traces and the extraction of the DIS 
limit are discussed in the Appendices. The parameters 
of our wave functions for both solutions are given in Ta- 
ble [VT] with the functional form of the S-state described 
in Eq. ([12} and the P and D-states in Eq. ( |4.10[ ). 

We find two solutions that fit the functions f(x) and 
gi(x), but only one of these (solution 2) also gives a sat- 
isfactory description of 32(2;)- This solution has a large 
D-state probability of (0.59) 2 ~ 35% and a small, but 
important P-state probability of (0.08) 2 ~ 0.6%. The 
other solution, solution 1, unfavored by the existing data 
for 32(2;), has a P-state probability of (0.43) 2 ~18% and 
a D-state probability of (0.18) 2 ~3%, values more in line 
with previous expectations [H HOI [41] . 

The observation that the solution with a large D-state 
exists and fits all of the data is perhaps the most impor- 
tant conclusion of this paper. It opens the door to new 
possibilities, and invites a more serious study of the pos- 
sible dynamical origin of the D-state components of the 



nucleon wave function. Furthermore, the fit requires that 
the u and d quark D-state components have a different 
sign as well as a different shape. 

The role that the data for (72 play in this conclusion 
focuses once again on the importance of accurate mea- 
surements of this small quantity. If the data on g 2 were 
to become sufficiently robust to allow the separation into 
u and d quark contributions, it could, from our point of 
view, be used to decisively fix the size of the P and D 
components of the wave function of the nucleon. 

Of course, many issues remain. First, we emphasize 
that our conclusions are sensitive to both the proton and 
the neutron data, but the extraction of neutron data is 
not without errors any uncertainties. And, the modern 
fits to both the unpolarized and the polarized data rely 
on QCD evolution equations to evolve the data to Q 2 > 1 
GeV 2 from a fitted phenomenological distribution defined 
at Ql = 1 GeV 2 . This suggests that Qq = 1 GeV 2 is a 
good place to match a quark model calculation to QCD, 
as we have done. However, the details of our fits would 
change if we chose to fit the quark model to QCD at 
a different Qq. Some investigators chose to match the 
distributions at much smaller Qq, even as low as Q\ = 
0.16 GeV 2 02]. The best choice of Q\ and even the 
validity of the procedure itself, is, in our view, still an 
open question. 

The orbital effects which we calculate refer to effective 
contributions which also include, inevitably, some contri- 
butions from glue. Thus our results may not be incon- 
sistent with recent lattice QCD calculations that suggest 
large contributions from the glue |43j . This cannot be 
discussed further until our model is used to compute the 
energy-momentum tensor and study the spin sum rule 
(1.2). This topic, alluded to very briefly in Sec. II G 



requires future study. 

We have not done a systematic fit to the data by min- 
imizing x 2 ■ If we did so, the fits would certainly im- 
prove, and we would be in a better position to assess the 
uniqueness of our final parameters. Until this is done, the 
parameters given in Table |Vl| must be considered prelim- 
inary. 

How realistic is such a large D-state probability? First, 
it must be emphasized that, because we did not fit all the 
wave function parameters including np and no in a sin- 
gle search, it is very possible that a good solution with a 
smaller D-state exists (and we think this is quite likely). 
So, without a systematic x 2 fit, we do not know precisely 
how large a D-state probability is required by the data. 
Here we present evidence of a large D-state. For compar- 
ison, note that the D-state probability in 3 H varies from 
about 7-9%, depending on the NN force model used (see 
Table V in Ref . 44 , for example) and at least four times 
smaller than our estimate. (The P state probability in 
3 H is also quire small, less than 0.16%; if this were also 
increased by the same factor of four as the D-state, it 
would compare with our result for the nucleon.) In the 
three nucleon bound state the D-state arises primarily 
from the large tensor force which is a feature of the one 
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pion exchange (OPE) part of the NN potential. If the 
force between two quarks is described by a combination 
of one gluon exchange (OGE) and a confining interaction 
(and/or OPE, as some models suggest), and the confining 
interaction is a mixture of scalar and vector components, 
then a tensor force can arise from both the OGE and 
the vector part of the confining interaction |45j . It may 
be that this mechanism can produce a large D-state, but 
Isgur, Karl, and Koniuk ,46J, using the harmonic oscilla- 
tor model of Isgur and Karl |47) find a very much smaller 
value. On the other hand, our findings are in line with the 
results of Refs. [JH1 US] on the importance of higher or- 
bital angular momentum components in the baryon spec- 
trum, affecting the nucleon mass by about 100 MeV, and 
producing inversion of the relative positions of positive 
and negative parity nucleonic excitations (unfortunately 
the breaking down of the contributions from each of the 
different partial waves included in those calculations was 
not reported). All that we can say here is that a large 
D-state can explain the DIS observations. A dynamical 
calculation is needed to determine whether or not a large 
D-state is possible. 

If the correct explanation of the DIS data is closer to 
solution 1 with the expected P-state and a much smaller 
D-state, then spin puzzle is easily explained by a nucleon 
wave function with modest angular momentum compo- 
nents. Otherwise, the spin puzzle leads directly to an- 
other puzzle - what are the dynamical origins of such a 
large D-state, and what effect does this have on other 
nucleon observables? 



the operation C, where 



7m' = 7m 



(Al) 



Since this operation does not depend on the quark or di- 
quark variables, it may be applied after the cross section 
has been calculated, allowing the construction of the full 
(conserved) current from the unsubtracted current 7^. 

For example, if the cross section arising from the un- 
subtracted current is written in terms of a general oper- 
ator O 



(A2) 



then the full result is obtained by applying the operation 
(lAlj) to both sides of Eq. dA2l). On the l.h.s. we get 



7m 



<1 



O 



(A3) 



which is precisely the result for the full current (2.4 1. On 
the r.h.s. we get 



^M v M v rvv iLV 



qM WW qu{q v w llu .) 



WvW q v WW) 



q 



(A4) 



which is conserved. Furthermore, if any terms propor- 
tional to either or q v arise from the unsubtracted cal- 
culation, they my be ignored because 
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Appendix A: Calculations of the structure functions 



1. Simplifications 



The full quark current is (2.4), but, because the sub- 



traction term ^q^/q 2 can be reconstructed from the first 
term, it is sufficient to calculate the 7 M contribution (re- 
ferred to as the unsubtracted contribution) only. To prove 
this, note that the full current can be constructed from 



C/ V - 0, (A5) 

1. e. they vanish once the full result is constructed. 

2. Contributions from squares of S, P, and D-states 

a. Isoscalar diquark contribution 

The only contribution from the D-waves to the 
isoscalar diquark term comes from the square of the i[> D '° . 
Together with the S and P-state terms this gives 



(e 2 q )° 
8M 



W, 



1=0 



ss 
11 



p'k v 



(A6) 



where m q + ft' ~ ft' , a factor of 1 /2 from the wave func- 
tions has been moved to the factor multiplying the double 
integral, a factor of \k 2 \ — > —k 2 coming from the defini- 
tion of ^i D ' is included in I DD , and the integral is 



p'k 



d 3 p'd 3 k 
(2^) 2 4e q E s 



S 4 ( P ' + k-P-q) 
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p! S = tr 



I pp = tr 

1 vp, — lL 



A^S + (ml-p' 2 )8 + (m 2 a -k 2 ). (A7) 
The S and P-state traces are 

K7i> f'l/J. 

where we use the shorthand notation 

V s = (M + f)(l + 1 5 $). (A9) 

The operator X DD is obtained by summing over the po- 
larization vectors e and £ using the identity (which holds 
for all polarization sums defined in the fixed-axis repre- 
sentation [24| ) 

Aq /3 = ^2( £ ^)a( £ *A)l3 = -9<*P ■ ( Al °) 
A 

Then, using 

n 

= -±%%, -^Dp, (P,k), (Aff) 
we obtain (including the factor of — k 2 referred to above) 

C = i fc 2 E( £ A) Q '( £ A) Q G Q ^'(fc,c:)G^(fc,c™) 



|tr[7^7 5 7,/ 7/i 7 5 7 ^Ps 



- -b A T g - J- k 2 T D 

vfi 20 v[i 



where we introduce two new standard traces 

la lv f'ln I a Vs 



J 9 = -tr 



T u = D /tr 



I 01 Ivf'l^l 01 Vs 



(A12) 



(A13) 



The trace P^ can be expresed in terms of the others. 
Using k a i a = K, 



pp 



I" = I 

V{1 ^f/J, 



(A14) 



Only the terms Vs -> f + Mj 5 $ contribute to all of 
these traces, and I pp and I 9 can be simplified by moving 
one of the I (or 7 a ) factors through the operator Vs, 
giving 



jpp 



= -k 2 tr 



Ivf lix-i 



-2M (S ■ k) tr 



fc 2 Mtr 



Ivf 7^7 



Ivf 7^7 



-k 2 1^ - 8M(S ■ k) ie^ aP fc> ,/3 , (A15) 



I 9 



-tr 



Ivf 7/J 



Mtr 



Ivf 7^7 



-§Mtr 



Ivf 7/x7 



.1 tSS _ 2 rSS 
3 vfi 3 i±v ! 



(A16) 



where the reader should be careful to distinguish 7^f 
from P™. Putting this all together gives 



tDD 



J- h 2 T PP _|_ 13 1-4 
20 1 vfi ~ 60 a 



.S.S" 



j-2 r 



(A17) 



These considerations show that all of the traces depend 
on only two elementary traces, I ss and J pp , with 



/^f = I{p', P) - 4Mie MVa p p' a S 13 (A18) 



where 



I( P ',P)=tT 



Ivf 7m- 



= 4^P„ + - <W ' P')]- (A19) 
Because q — > oo, we may replace p' — > q in the last terms 



of (A18) and (A19|, but not i n th e terms proportional 
to p'^ or p' v because, through ( A5 1 , the q contributions 
to these terms will vanish. To reduce these we use the 
identity 

i J dfl k± p'„ = A X P^ + A 2 q^ -> AiP„ + q M (A20) 



where the nonleading term, A\ , given in Eq. ( B5 ) , is 



needed only when the leading q^ contribution vanishes. 
To evaluate the last term in ( A15 ), first note that p' = 



q + BP-k where B = [1 - (k ■ P)/M 2 ]. Hence the p' may 
be replaced by q + BP and we may use identity ( B6 ) to 
reduce it. In doing the reduction, note that 

-8M{S-k)ie^ aP k a {q + BPf 
-> -8MC X ie^ a/3 S a (q + BPf 

-8MC 2 (S ■ q) ie^ a0 (q - ~P) a (q + BPf 
-> 8M 2 d T x - 8 Af 3 v C* 2 T 2 ■ ( A2 1 ) 

With these simplifications, the two standard traces be- 
come 



jSS 4 



I pp =4 



2M 2 Ai 



P M P, 
M 2 



(P ■ p')9w - m 2 i, 



1k 2 M 2 A 1 ^)+^(P-y) V 
— M 2 (fc 2 - 2 Ci )Xi - 2 A/ 3 f C 2 X 2 j . (A22) 



Using the relations (valid in the DIS limit and derived 
in Appendix |b|) 

k 2 - -k 2 
(P ■ p') -> Mv 
A\ -> x 

C x = -ifci = -ik 2 (l-z 2 ) 
fc 2 -2d = -k 2 z =-|k 2 [l + 2P 2 (z)] 

v 2 G 2 = c 2 -> k 2 (|z 2 - i) = k 2 P 2 (z) (A23) 



21 



where k — {0, k} and z — cos 6, the three traces from Eq. 
(A6) reduce to 



jSS _ 



t pp 



4M 
4Afk 2 



2M.r [ P ^ ] -vg^-MIx 



M 2 



, , 2M 2 
-P 2 {z){ \MX X X 2 



M 2 



%MX X 1 2 

6 v 



7^=4Mk 4 
These traces are all of the general form 

z l9fuy + z 2 



(A24) 



(A25) 



Including the factor of l/(27r) from (2.1), the invariants 



then have the general form (with / = or 1) 

M n 

i p'k 



Zl n 2 L ^ 2 L (P,k) (A26) 



where Z t = {Wi,W 2 , Gi, G 2 } (with i = {1,4}) and z 2 is 
the coefficient from (A25). 



Combining all of these results gives 



2MxWl=° = vW{= Q 



G[=° = (el) 



,2Mvx 

4-7T 

M 



4?r 



p' k 

p'k 



2 /2 , 2 1 2 / 2 i 2 1 4 / 2 

n s tpg + n P k tpp + n D k ip D 



2/2 1 2 1 2 / 2 1 2 1 4 / 2 



-P 2 (z)[ 2 3 n 2 P k 2 i: 2 P + ^n 2 D k 4 ^ 



M 



G 



1=0 



<e 2 q ) 



2tt 



3 '"P " f P i 30 



where ipL — iPl(P> k), and 



Gi = Gi + ~j^?G 2 . 



(A27) 



(A28) 



Hence, the isospin zero structure function G\ becomes 



4?r 



p' k 



2/2 1 2 i»2 / 2 1 2 1 4 / 2 

rip k ip p — i%n D K ipi 



+P 2 (z) | np k^p + j s n 2 D k>£ • (A29) 



Note that the W\ and structure functions satisfy 
the Callen-Gross relation, which will hold for all contri- 
butions calculated below. 



b. Isovector diquark contribution 

Next, look at the contributions from the 7=1 diquark. 
Here there are four contributions: from the squares of 
^ r S,i j ^p,X j ^rAl ) anc i All of these involve sum s 

over the diquark polarization A evaluated using 
Carrying out the spin sum and removing the 7 s 's (which 
introduces another minus sign) gives four traces 



(e 2 ) 1 

w i=i = \W_ 
^ 8M 



T PP 



{n 2 ^ 2 (P,k)J^ + n 2 p^ 2 P J^ 
p'k K 

+n 2 D i/>% (P, k) (j# + J 2 - 2 ) } (A30) 
where, using ^f a f + f^j a =0 and 



§ k 2 g aa ' + \k a k a ' 



(A31) 



these traces are 



tSS _ 



-ftr 



J, 



pp 



3 vp- 

-itr 



la lu / 7m 
1 tSS , 2 rSS 

"" "T" 3 1 pu 



= -I 9 





IjPP , 2 
3 I'M ' 3 


Tl,l 


5 K i ^ 


/ 





2 .2 jPP _ 2 ,4 rg 
5 5 * J iyp ' 



(A32) 



In the previous section, all of these tr aces w ere expressed 
in terms of the two standard traces (A22). Adding the 



two D-state contributions together, and simplifying, gives 



t dd 



2 ,2 jPP _3 k 4rg 



2 p jPP 1 1 1.4 



AMk 4 i 2Mx 



P^^Ml.-^Ll,} 



vg^v 1 3 
2Af 2 



\Mli 



(A33) 



Inspection of the results for J ss and J pp shows that 
contributions to the unpolarized structure functions com- 
ing from isovector diquarks are equal to the isoscalar 
ones, while for the polarized structure functions, J^/f ~~ ^ 
— jjl„u and — > Using these arguments, we 

obtain the following results 



2MxW( =1 = vWi =1 



{elf 



Mvx 
2tt 



p'k 



2 / 2 1 2i2/2 1 2 1 4 / 2 

n s ip s + n p k tpp + n D k ip D 
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1=1 



{ q} 12tt 



p'k 



n 2 s ^ 2 s -\n 2 P Y 2 ^ 2 P +n 2 D Y^ 2 D 



InlkVp-fr^kVI, 



(A34) 



Extracting Gi gives 

Gl " ^ 127T 



p' k 



+8P 2 {z) I n 2 P kVp - W D k 4 ^ 



(A35) 



Now add the isospin and isospin 1 parts, separating 
them into the u and d quark contributions (for the pro- 
ton) following the discussion which lead to Eq. (2.26|. 
Collecting terms gives 



2MxWf = vWl = x 



Mv 



M 

6-7T 



p' k 



2e 2 u N u + e 2 d M d 



p'k 



4e 2 (^) 2 -e 2 (^) 2 



§e 2 k 2 (^) 2 -I e 2 k 2 (^) 2 



2e 2 k 4 (^) 2 + e 2 k 4 (^) 2 



+n P P 2 (z)[f e 2 k 2 (^) 2 -|e 2 k 2 (^) 2 
+^^W[ie 2 k 4 (^) 2 + | e 2 k 4 (^) 2 



M &2 



M 



P2(Z) 



p'k 



|e 2 k 2 (^) 2 -I e 2 k 2 (^) 2 _ 
ie 2 k 4 (^) 2 + | e 2 k 4 (V^) 



(A36) 



where 

N q ee n 2 (iP s q ) 2 + n 2 P k 2 «) 2 + n 2 D k 4 «) 2 (A37) 

is the density factor that appears in the unpolarized 
structure functions. 

3. Extraction of the DIS limit 

Before calculating the interference terms, we extract 
the DIS limit of the above results. The expressions for 
all of the structure functions are covariant, but it is con- 
venient to evaluate them in the laboratory system, where 
the two external four- vectors are 

P = {M, 0, 0, 0} 



2Mx 



,0,0, \ Q 



Q 4 



AM 2 x 2 



(A38) 



where x is the usual Bjorken scaling variable. For the 
integration of functions of |k| and k z = \\c\z, we can write 
in the DIS limit 



p'k 



d 4 k 

d 3 k 



S + (m 2 q -p> 2 )5 + (m 2 -k 2 ) 



Mx 



k 2 dlkl 



dzS([M(l-x)- E s + \k\z}). 



Scaling all momenta by the nucleon mass (so that |k| 
Mk) the double integral becomes 



p'k 



M 2 x ndn 
~W Jo 



4tt£L 



ndz S ((1 — x) — E K + kz) 



M z x 

~Q 

rris 
v 



ndn 
A^E~ 

dx_ 
16tt' 



dz S (zq — z) 



(A39) 



where \ was defined in Eq. ( 2.35 1, E K = \Jr 2 + k 2 , z = zq 
with zq defined in (2.391, and the lower limit of the % 



integral (which occurs at | zo I = 1) wa s given in Eq. (2.411 
To obtain the final expression we use 



16tt 2 



dx = 



Mv 



C 



p'k 



(A40) 



When inserted into ( A26 1 , the structure functions have 
the form 



i 2I Mm s 
1 { q) 2v{An) 2 



d X z in li>l{P,k) (A41) 



C 



where the Zi are the factors that appear in Eq. (A25l 
For the proton this is 



vWii = 2MxW[ 



v _ 



2e 2 u f u {x) + e 2 f> d {x) (A42) 



where [after the SP interference terms have been calcu- 
lated, this /' is replaced by the f q given in Eq. (2.31)] 



K = n 2 s f q S + n 2 P fP + n 2 D f q D 



(A43) 



with the individual quark distribution functions f q (x) 
given by 



Mm s 
(4tt) 2 



(A44) 



where now k = |k|, and L = {S,P,D} with k 21 = l,k 2 , 
or k A for S, P, or D-states, respectively. Up to the normal- 
ization factor, treated differently in this paper, the struc- 
ture functions f q (x) were already obtained in Ref. I15|. 



The interpretation of Eq. (jA44|) was discussed in Sec. 
above. 
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The other proton structure functions become 
vGl = gl{x) = n 2 s U4elfi(x) - e 2 f d {x)] 



,2 1 



-n 2 D \[2elf°{x) + elf°{x)] 
+n%U ie l9n(x)-e 2 d g^(x)] 
+n 2 D ±[29e 2 u gZ(x) + 16e 2 d g d D (x)] 
gl(x) = -n 2 P \[Ae 2 u g? (x) - e 2 d g%(x)\ 

-«dw [29 e 2 u g^(x) + 16 e 2 g°{x)\ ( A45) 



where the new structure functions, both of which depend 
on P2(zq), are 



9q(x) 

9a {X) 



(4tt) 2 
Mm s 



d X P 2 (z )k 2 [^( X )} 2 
d X P 2 (z )k^( X )} 2 , 



(A46) 



c 



where zq was defined in Eq. (2.391 



4. DIS limit in light cone variables 



In light-cone variables the rotational symmetry is bro- 
ken so the variable X , which previously depended only on 
k, now depends on both x and k\ = Mm s i] 



X{v,x) 



k+ + k. 



- 2 



x 



m s 
1-x 



-2 = C 



1-x 



(A51) 



However, the integral (A50) can be immediately trans- 



formed into ( A39 ) , showing that the representations 



are equivalent. We used the light cone form (A50) in 



Ref. [15] , but prefer ( A39 ) for this paper because of the 
interpretation of its meaning, as discussed in Sec. |II E| 
above. 

In order to finish the comparison, note that, in light 
cone variables, the weight functions that appear in (2.28) 
are 



fc_) = 2 m 6 



kz Q = k z = \{k + 

k 2 (l-z 2 ) = kl=Mm s y(x-C) 
k 2 P 2 (z ) = k 2 z -\kl, 

where, as elsewhere, y — 1 — x. 



Xl 



(A52) 



The DIS limit can also be evaluated in light-cone co- 
ordinates. In our notation, an arbitrary four-vector v is 
written 

v = {v + ,v_, vjj v ± =v ±v 3 (A47) 

so that the scalar product is 

v ■ u = v^u^ = \(v + u- + «-it+) — vj_ • uj_ • (A48) 

In the DIS limit, the four-vectors P and q, in light-cone 
notation in the laboratory frame are 



P = {M, M, 0} 



(A49) 



Using this notation the double integral can be reduced 
to an integral over k\. Defining fc_ = M(l — y), the 
integral becomes 



p'k 



d A k 
(2^ 



5 + {m 2 q -p' 2 )5 + {m 2 s -k 2 ) 



J^\J dk^J dk+S(m 2 s + kl-k+k^ 



4tt 2 _ 



dk]_ 
~8V 



dy 



(i-w) 

dk\ 



Mv J 16tt(1 - x)' 
where the S(y — x) fixes fc_ = M (1 — x). 



(A50) 



Interference Terms 



SP interference 



The contribution to the SP interference term from 
isoscalar diquarks is 



(eg) 
8M 



W, 



7=0 



n P n s ^s(P,k)yj P (P,k)I^ (A53) 



p'k 



where the new trace is 



SP 



tr 



(A54) 



and we have been careful to separately write the two 
contributions from the overlap of the final P-state (the 
first term) and the initial P-state (the second term). As 
it turns out, these terms are not identical. To reduce the 
calculation, note that the trace now picks up the terms 
with an even number of 7 matrices in Vs 



V s -> M + f>-y 5 $. 
Moving the k through Vs reduces the trace to 



(A55) 



1% = 2Mtr 



7p 



2{k ■ 5)tr 



8M[K P '^+p% - (k-p')g^ 
+8i{k-S) e^ apP ' a pP. 



(A56) 



In order to expand this into the four independent DIS 
structure functions, we recall that p' = q + BP — k (with 
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MB = M — E s ), and average over the directions of kj^ 
using the identities (B2| and (B6|. We get 



= 8M [.Bi (q v (q + BP),, + (q + BP) v q^ ' 

-ICxg^v - 2C%q li q ll - (k ■ p')g^ 
+8M 3 {B l ~C 2 )l 2 + 8i Ci e^ aP P a S f3 (A57) 

Dropping all terms proportional to to q^ and q v , the un- 
polarized terms reduce to 



rSP 



-8M g^[2d + (k ■ p')] 
P,P 



-8M 



2MvB x B - 2(d - v z C 2 ) (A58) 



and in the DIS limit the coefficients can be simplified: 

2C*i + (k ■ p') -k\ - k z q z + k 2 ->• -k z v 
2MvBB 1 - 2(Ci - C 2 v 2 ) -> 2(M - E s )k z + 2k 2 z 

2k z (k z +M- E s ) = 2k z Mx . (A59) 

where the light-cone relation fc_ = E s — k z = M(l — x) 
has been used. Substituting, gives 



SP 

V ft 



=8M 



k z vg tlv -2k z M: 



AP 



(A60) 



The polarized terms contain an invariant not in the 
canonical form, ft can be re-expressed using the identity 



(C3|, giving 



jSP 



■MZClEy) 



+T 2 [8M 3 {B l -C 2 ) + 8C l E 2 ]. 

In the DIS limit these coefficients become 
k 2 

8Ci£i ->• 8M -i 
Ax 



(A61) 



8M 3 (B 1 - C 2 ) + 8dE 2 
'k z 1 



8M i 



k 2 -±k 2 



AM 3 % 



8M 



2 r 



Mk 7 



k\ 
Ax 



giving 



[S P 



= 8AP 



AMx 



M 



k z - 



AMx ) 



(A62) 



(A63) 



Using ( A25 ) and ( A41 1 , the isospin zero contribution 



to the SP interference term is then 

vW 2 I=Q = 2MxW(=° = -n P n s {e 2 q )° x h° q {x) 
vG[- =-n P nsl{e 2 q fh° q {x) 



M 



Gi=°= -n P n s \ [e\f h\{x) - h° q (x) , (A64) J„„ 



where we introduce two new structure functions involving 
the overlap of the S and P-states 



h° q (x)= 
hl(x)= 



Mm a 



(4tt) 2 J c " v : ' ! x < 



Mm s 



C 



(A65) 



It is now an easy matter to compute the isovector con- 
tribution to the SP interference term 



W. 



i=i 



(eg) 1 
8M 



n s n P xP s (P, fc)Vp(P, k) J% (A66) 



p' k 



where, removing the 7 s 's and summing over the diquark 
polarization, gives 

■C = |tr [(7a ^/7„ 7 Q + % T)Vs] ■ (A67) 

Using the identity 

%{M + f-i 5 $W = 3M - ff$ (A68) 

we see immediately that the unpolarized part of J SP 
equals the unpolarized part of I s p , while the polarized 
part of J SP is — J of the polarized part of I s p . This can 
be be summarized by the relation 



J, 



SP 



1 tSP , 2 rSP 

3 1 up ' 3 i ^ii/ 



(A69) 



previously encountered for J bb and J in Eq. (A32| 



Adding the two isospin contributions, and separating 
the u and d quarks, gives the following result for the 
proton structure functions 



vW%{x)= 2MxWf 



g{(x)= -npns[i e t h °u( x ) ~l e d h di x )] 



2n P n s x[2ei h»{x) 

1 

3 

g%(x)= -n P n s (j e 2 u [h} u {x) - h° u {x)} 
[h d (x)-h° d (x)} 



1 P 2r 

3 e d 



Mix)] 



(A70) 



These are combined with the S and D-state contributions 



as reported in Sec. IIP 



b. SD interference 

Only the fy D ' 2 D-state component can interfere with 
the S-state components, and this involves an isovector 
diquark, giving 



W, 



1=1 



(e 2 ) 1 
8M 



p' k 



a S DMP,k)i> D {P,k)J^ (A71) 



where qsd, defined in Eq. (2.32), includes some factors 



from the D-state wave function. Summing over the po- 
larizations of the diquark and removing the 7 s 's, gives 
the new trace 



SD 



-§zr a tr 



1^ (A72) 
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where the two identical terms (with the D-state in ci- 
ther the initial or final state) have been c ombin ed, and 
the trace I D was encountered before, Eq. (A13). Using 
(|A14|, ([AlBl and (IAIdI), 



t d _ 11.2 

1 VH ~ 3 



rSS _ tSS 



8M(S ■ k) ie^ a p k a p'P 



{ci-fk 2 )- 



= 8M 

8M 2 k 2 P 2 (z) 
This contributes 

~il=0 



Ci - ^k z )1 1 - MvC 2 Ii 



M. 



(A73) 



M 



G 



\a D {e 2 q ) l d q {x), 



(A74) 



where the new structure function is 



d q (x)= 



Mm* 



(4*0 ^ c 

The proton structure functions are 



d X P2(z )k 2 ^( X )^ q J (x)- (A75) 



2 

qCL D 



e 2 u d u {x) + 2d d (x) 



g%(x) = \a D 



e 2 u d u (x) +2d d (x) 



(A76) 



c. PD interference 

As for the SD interference, only the ty D ' 2 term will 
interfere with the P-state, and the trace is 



(e 2 Y 

1=1 _ \ C q> 



w, 



8M 



a PD 4>d(p, Wp(P, k) J™ (A77) 



p' k 



from the D-state wave function. Summing over the spin 
of the diquark and removing the 7 5 's gives 



jPD _ 1 niffl't 



1L2 



k z tr 



(la lv f 7m ^ 7a' + 7a ^ li> f 7/i 7a' ) Ps 

{% j v f' 7^ + 7„ f' 7 M iCj Vs 



where I SP \ P was given above, Eq. (A63|. This term con- 
tributes to both G\ and G 2l giving 



M 



G 



i=i. 



apD(e 2 q y h*(x) -h 2 (x) 



(A81) 



where the result has been expressed in terms of the new 
structure functions 



h\{x) 



Mm* 



(4tt) 2 J ( ' " ' •' 



dxk 3 z i fl (x)t P (x) 



K(x) 



Mm s 
(4tt) 2 



k 4 (l 



4Mx 



<(X)<(X)- (A82) 



Evaluating this for the proton gives 

gl(x)=la PD [e 2 u h 2 u {x) + 2e 2 h 2 d {x)] 
g% (x)= § a PD (e 2 u [hl (x) - h 2 u (x)} 

+2e 2 [h 3 d (x)-h 2 (x)] 



(A83) 



Note that the structure of thess terms is similar to the 
SP interference terms, except for the charge weightings, 
which parallel the SD interference terms. 



Appendix B: Details in the derivation of the DIS 
formulae 

First, evaluate the average of k a over the directions of 
kj^. Using the fact that the wave functions are indepen- 
dent of the direction of k^, and that P ■ k — 0, we use 
the fact that the integral can depend only on 



9 



(P ■ q)P 

M 2 



vP 



where a PD , defined in Eq. (2.32), includes some factors which gives 



1 
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dfl k± k a = B 1 q a . 
The coefficient B\ is therefore 



B 1 = 



k-q _ M(k ■ q) - v{k ■ P) 



q 2 



Mq 2 



|q| 



k' z 
V 



(Bl) 



(B2) 



(B3) 



(A78) 



where, in this term, Vs — > M + fi b $. The unpolarized 
terms cancel, and using the identities 

% ff$ + = 2{k • SlfV 

ll a fl b $l a + l a fl 5 $T I = -2(k ■ S)fj 5 (A79) 

this trace reduces to 



where the last experssion is the result in the DIS limit. 
The average over p' a = q + P — k = q + BP — k [where 
B = 1 — (k ■ P)/M 2 ] follows immediately 



2tt 



BP a - B x q a 



A x P a + A 2 q 



(B4) 



where A± and A 2 are 



J^° = fk 2 (fc-S)trk/ 7 ^7 



4 k 2 tSP 

9 vii 



(A80) 



A 2 =l-B 1 



1 



(B5) 
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The second integral we encounter is the average 
1 



2tt 



dft k± k a k/3 = Cig a p + C 2 q a qp 



(B6) 



where g was defined in Eq. (2.17). The simple form (B6) 



follows from the conditions that the contraction of P into 
either index must give zero. The coefficients are found 
from the relations 



k 2 = 3Ci + C 2 f 
(q ■ k) 2 = C,q 2 + C 2 q 4 



(B7) 



giving 



C 2 



q 2 k 2 - (q ■ kf 
2q 2 

3(g ■ k) 2 - q 2 k 2 
2q* 



(B8) 



Noting that q- — — q and k = — k. we can write 



Ci = -\k\ = -ik 2 sin 2 i 



C 2 = 



1 



3Z,2 _ 11,2 
2 A, Z 2 



In the DIS limit 

<f = -(q 2 

we may write 

C 2 -> ■ 



Af 2 



p _ 1l2 



-> -1/ 



1 



i- 2 _ It 2 



To express C 2 in terms of fcj_, we note that 
fc_ = M(l - x) 



= E s -k z = 
Solving for k z gives 

k 



m 2 + k 2 . 



2M(1 - x) 



m 2 + fc 2 ^ + fc 2 — fc 2 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 



and hence 



C 2 



(m 2 + fcj[) 2 
4M 2 (l-.x) 2 



i(m 2 + 2fci) + iM 2 (l- a; ) i 



■ c 2 . 



(B14) 



Even though C 2 is very small, it cannot be neglected 
because it multiplies a large term. 
Note that 

k 2 - 3Ci = C 2 q 2 -> -c 2 = -fc 2 + |fc 2 . (B15) 

Appendix C: Identity for the reduction of the 
hadronic tensor 

Here we prove an identity needed for the calculation of 
the P-state contributions to the hadronic tensor. This is 



= E l l l +E 2 l 2 . (CI) 



To prove this identity note that both sides conserve 
current and are antisymmetric; therefore the coefficients 
Ei can be determined by requiring the the projections 
P^ T„ v and S» Tpv be equal. This gives 



Ey- 



E 2 = 



M 

M 3 _ M 2 
~Q 2 ~ 2^? 



(C2) 



In our calculation we were able to ignore terms pro- 
portional to q^ or q v . If these terms are dropped, the 
identity becomes 



ieu»«8P a S f} ^ E 1 l l + E 2 l 2 



(C3) 
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